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Abstract 

We obtain limit theorems for the row extrema of a triangular array of zero-modified geometric 
random variables. Some of this is used to obtain limit theorems for the maximum family size within 
a generation of a simple branching process with varying geometric offspring laws. 
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1 Introduction 

It is well known (Anderson (1970)) that the geometric law is not attracted to any max-stable law and 
hence maxima of independent geometric variables cannot be approximated by a max-stable law. Con- 
sidering triangular arrays of zero-modified geometric laws allows adjustment of the zero-class probability 
independently of the success probability parameter, thus opening the possibility of approximating row 
maxima and minima by simple explicit laws. Similar results for Poisson laws are given by Kolchin et 
al. (1978, §2.6). Our motivation is closer to that of Anderson et al. (1997) who study the Poisson and 
other laws. They exploit the normal approximation to the Poisson law with large mean to show that row 
maxima are approximated by the Gumbel law under certain conditions. We obtain corresponding results 
which emanate from the exponential approximation to the geometric law when its mean is large. 

For each n = 1,2,... we let v n be a positive integer and {X,(n) : i — 1, . . . , v n } be independent 
random variables with the same zero-modified geometric law 

P{X l (n) = {)) = l-g n fc P(X l {n)= ] )=a nPn {l- Pn y-\ (j = 1,2,...), (1-1) 
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where < a n < 1 and < p n < 1. The mean for row n is a n /p n and the distribution function is 



where [x] denotes the integer part of x. The standard geometric law corresponds to a n = 1 — p. 
In the next section we prove limit theorems as v n — > oo for the row extrema and range 



and we give examples showing our hypotheses can be satisfied. More specifically, we find conditions 
which ensure the row extrema converge in probability to infinity, and show in Theorems 1 and 3 under 
a further condition that normalized versions have non-defective limit laws. Theorem 5 demonstrates 
their joint weak convergence, and Theorem 6 shows that the range is asymptotically proportional to the 
maximum. Lemmas 1-4 provide context by exhibiting possible behaviours of the extrema under differing 
assumptions. 

Our results for row maxima are used in §4 to obtain corresponding limit theorems for the maximum 
family size, again denoted M n , in the n-th generation of the simple branching process where the offspring 
law for individuals in generation n — 1 is the geometric law (1.1). Of course this is precisely the case of 
a varying fractional linear offspring law which has previously been studied by Agresti (1975), Kciding 
and Nielsen (1975), and Fujimagari (1980). Maxima of random variables defined on the classical Galton- 
Watson tree have been studied by Arnold and Villasenor (1996), Pakes (1998), and Rahimov and Yancv 
(1997,9). We show in Theorems 9-11 that results from §2 transfer to the branching process setting through 
conditional limit theorems (Theorem 7 given in §3) for the generation sizes Z n . These latter results seem 
to be new, and they are the strongest possible assertions, which can be made within our restricted class 
of offspring laws. 

2 Behaviour of row extrema 

To set our main result for A4 n in context we begin with some elementary descriptions of its behaviour. 
Observe that the distribution function of M n is H n (x) :— F£ n (x). 

Lemma 1 Let lim n _>oo v n = oo. (i) M n if u„a n — ► 0. 

(ii) If ^i/„o„ < oo then P(M n > i.o.) = 0, and if the rows are independent and ^v n a n — oo then 
P(M n > i.o.) = 1. 

(iii) Let < e < 1. If the rows are independent and u n a n > logn + (1 + e) log(logn) for all large n then 
P{M n = i.o.) = 1. If v n a n < logn + log(logn) then P(M n = i.o.) = 0. 

Proof For (i) observe that H n (x) — > 1 for all x > 0. The remaining assertions follow from the Borel- 
Cantelli lemma and elementary estimates of P(A4 n > 0) < 1 — (1 — a„)"") and P(A4 n = 0). O 
The quantity a n = log(^„a„) is important to our further considerations. The following limit theorem 
is easily proved. 




(1.2) 



Ai n = max Xi(n), 



Ki<v, 



= min Xi(n), & 1Z n = A4 n — 
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Lemma 2 Let linin^oo v n = oo. (i) If 

lim a n — a (—00 < a < 00) (2.1) 

and 



l im ^ W = ( eXP( " eQ(1 " P)H) if ^ ' (2-3) 
.0 if x < 0. 



lim p n =P (2.2) 

then 

1 " I, 
<ii) it 

lim a„ = 00 & lim sup p n < 1, (2.4) 

then A4„ 00. 

The limiting distribution function in (2.3) is non-defective if p > 0, and it is defective if p = 0. In the 
latter we have M n — > A^oc where 

P(Xoo = 0) = 1 - P(Xoc = 00) = G(-a), 

and G(x) = exp(— e _a: ) (—00 < x < 00) is the distribution function of the standard Gumbel law. 
Theorem 1 characterizes the rate of divergence to infinity under some further conditions. 

Theorem 1 Let rim n _,oo v n = 00. Assume that for some real c 

lim p n = & lim a n p n — 2c. (2-5) 

n — >oo n — >oo 

(i) If lim^oo a n = 00, then c > and 

PnM n - a n -i- A - c, 

where A has a standard Gumbel law. 

(ii) If lim„^oo a n = a (—00 < a < 00), then 

PnM n -4 (A + a)+. 

Proof Define 

X n = (x + a n )/p n & T] n (x) = V n a n {\ - p n ) [Xn] (-00 < X < oo). 

By using the expansion [x n ] = (x + a n )/p n — 5 n , where < 5 n < 1 is the fractional part of x n , it is easily 
seen that 

\ogrj n (x) = a n + [x n ] log(l-p n ) =a„+ r + a " {-p n - \p\ + 0{p 5 n )) = -x + 0(p n )- $a n p n (l + o(l)) . 

V Pn J 

The right-hand side converges iff (2.5) holds, and the limit is —x — c. Thus (2.5) is equivalent to 

lim 7] n (x) = e~ x ~° (-co < x < oo). (2.6) 



The distribution function G n (x) of p n M n — a„ equals [1 — u n 1 rj n (x)) l ' n if x n > 0, and equals zero 
otherwise. It follows that 

{G(x + c) if liminfn^oo x n > 0, 
if lirnsup^^ x n < 0. 

But x n — » oo for all real x iff a n — > oo, and assertion (i) follows. If (2.1) and (2.5) hold then c = and 
x„ — ► ±oo according as x > —a or x < —a, respectively. It follows from (2.7) that 

f G(x) if x > —a, 
lim G„{x) = < 

10 \fx<-a, 

and hence p n M n — ol max(— a, A), and (ii) follows. O 
Observe that the first member of (2.5) implies that if n 3> 1 then F n (x/p n ) w 1 — a n + a n (l — e~ x ), 

the exponential approximation mentioned in §1. Also note that Lemma 2(h) holds under (2.5). 

The assumptions of Theorem 1 can be realized. Let a n = e~ a v~ b (S > 0) and p n = (C > 0). 

Then (2.1) holds iff < 6 < 1, and then (2.5) holds with c = 0. The condition (2.4) holds if S = 1, 

and then we can admit any p„ — > 0. Now let S < 1 and choose p„ ~ A(log^„) _l » where A is a positive 

constant. Then (2.1) still holds and 

(1 - S) logv n 
a n p n = A- 




(logf„)< 

Then Theorem l(i) holds with c = if ( > 1 and with c = (1 - 5)^4/2 if C = 1. The case ( < 1 is an 
instance of (2.5) where c = oo. In this case Theorem l(i) suggests that the limit law is concentrated at 
— oo, and indeed this is true. In fact there is no affine transformation of A4 n which has a non-degenerate 
limit law. However we have the following large deviation estimate, 

log P(M n > Xn) = _1 
n^oo a n p n 2 

Further consideration of rj n (x), defined in the proof of Theorem 1, shows that if (2.1) holds and if (2.2) 
holds with < p < 1 then no affine transformation of A4 n has a non-defective limit law. The following 
result, generalizing the direct assertion of Theorem 2 in Anderson (1970) and with more explicit centering 
constants, shows that it is possible to stabilize the law of M n . The proof is similar to that for (2.6). 

Theorem 2. Set C n — — a„/log(l — p n ). If lim^oc a n = oo and (2.2) holds with < p < 1, then 

G(-y(x - 1)) = liminf P(M n - C n < x) < \imsupP(M n - C n < x) = G(-yx), (-oo < x < oo), 

n *oo n — >oo 

where 7 = — log(l — p). 

Parallel to Lemmas 1 and 2 we have the following results for the row minimum fi n , and they are easy 
consequences of its distribution function 

[0 if y < 0. 
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Lemma 3 Suppose that v n — » oo. (i) M n ^ iff <» - 0. 

(ii) If X^ a n" < 00 tncn P(m« > i-°0 = 0) an d if tnc rows are independent and J2 a n n — 00 then 
P(Mn > i.o.) = 1. 

(iii) If the rows are independent and ^2(1 — a^ n ) = oo then P(p n = i.o.) = 1, and if ^(1 — a^ n ) < oo 
then P(fi n = i.o.) = 0. 

Lemma 4 Let linin^oo v n = oo. Assume that 

lim z/„(l - a„) = /3 & lim v n p n = p. (2.8) 

(i) If /3 + p < oo, then 



lim K n (y) = 



1 - e-' 3 -*] if y > 0, 
if j/ < 0. 



(ii) If /3 = p = 0, then fi n oo. 

The limit law in Lemma 4(i) is non-defective if p > and we see that it defines a zero-modified 
geometric law. If p = then /i„ — ► /ioo where 

P( Moo = oo) = 1 - P( Mco = 0) = e^. 

Our principle result shows that if (2.8) holds with p = then p n can be centered and scaled to give a 
non-degenerate limit law. Set 

(3 n = -v n \oga n . 

Theorem 3 Let linin^oo v n = oo. Assume that for some be [0, oo) 

lim v n p n = & lim (3 n p n = 2b (2.9) 

and 

lim = G [0,oo), (2.10) 

then 

VnPnPn ^ (£ - P) + (2.11) 

where £ has a standard exponential law. 
Proof Let 

Vn = (y - l3n)/VnPn & Ipn (y ) = (a n (1 - p n ) [Vn] ) (-00<y<Oo). 

Observing that 

log^„(y) = -/?„ + v n (- — — - 5 n ) {-p„ - \p 2 n + 0{pl)) = -y+ \{3 n p n + v n 5 n p n (\ + o(l)), 
where <5„ is the fractional part of y n , it is clear that (2.9) is equivalent to 

lim ^ n {y) = e b -\ (-00 < y < 00). (2.12) 
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If (2.10) holds then 6 = and y n — > ±00 according as y > /3 or y < (3, respectively. Since K n (y n ) = 
1 — ip n (y) if y n > 0, = otherwise, it follows that v n p n [i n +/?„—> max(/3,£ ), and (2.11) follows. O 

The proof shows that (3 < 00 is a necessary condition for a non-defective limit law. The limit assertion 
(2.11) extends for lim^oo v n p n > in a manner similar to Theorem 2 as follows. 

Theorem 4 If lim^oo v n p n = £ € (0, 00) and (2.10) holds then 

1 - e~ y ~ ~ 6 < liminf P{v n p n \i n < y) < limsupP(f„p n/ u„ < y) < 1 - e~ y - f3 . 

n >oo n — >oq 

Observe again that (2.10) implies the second member of (2.9) with 6 = 0, and that a n — > 1. Hence 
our assumptions for (2.10) are precisely (2.8) with p = 0. In addition, p„logf„ — > if p = and 
a n = log^„ +o(l). Thus, the assumptions for Theorem l(i) are satisfied with c = and hence we have 

PnMn - logZ/„ -> A & V n p n [l n ->■(£- f3) + . 

Our next result extends this pair of weak limit statements to joint convergence, showing in particular 
that /j, n and M. n are asymptotically independent. 

Theorem 5 If lim^oo v n p n = and (2.10) hold, then 

(p n M n - logu n , v n p n fJ,„) -> (A, {£ - /?)+). 

Proof If x n = (x + log v n )/p n and y n = y/ v n p n , where x, y are real, then x n - y n = (y n log v n + v n x - 
y)/v n p n — ► 00. Consequently for any real x, y and n large enough we have 

A n (x,y): = P(M n < x n ,fj, n < y n ) = H n (x n ) - P(y n < fi n ,M n < x n ) 



Hn(Xn) 1 ' F n (x n ), 
The proofs of Theorems 1 and 3 show that H n (x n ) — » G(x), and hence f^^n) — ► 1, and v n F n (y n ) — > y+/3 
if y > 0, -> otherwise. It follows that A n (ar, y) G(x)(l - e^" ) if y > 0, -> if y < 0. O 

Our last result shows that the row ranges lZ n — M. n — fi n are determined by the row maxima. 
Theorem 6 If the conditions of Theorem 1 hold, then lZ n has the same limit behaviour as M n . 

Proof The assumptions imply that p n — ► 0. For any y > we have 

- log P(j> n [i n > y) = -P n - v n y + 0(v n p n ) -> -00, 
i.e., PnUn 0. The assertion follows from Slutsky's lemma. O 

3 The simple branching process 

Let (Z n : n > 0) denote the generation sizes of the simple branching process with varying geometric 
environments, 

Z n -1 

Z n =Y J *M) (« = 1,2,...) 
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where Zq — 1 and the -Xj(n) (i,n > 1) have the same geometric laws as in §1, and they are mutually 
independent. Thus X\(n) is a generic family size of a parent in generation n — 1, and its probability 
generating function (pgf) is 

(1 - g)jn + g 

/n(s) = (1 - .)r B + 1 (3 ' 1} 

where 
and 

m n = fn( 1 ) = a n/Pn 

is the mean n-th generation family size. 

As is well known (Harris (1963)) the pgf of Z n is obtained by functional composition, 

<P n (s) :=E(s z "\Z = l)= ( f>n-i(fn(s)). 

The group structure of Mobius transformations (3.1) permits the explicit determination 

(1 - s)A n + 8 
(1 - s)B n + 1 

where 

n n 
3=1 3 = 1 

M = 1 and 

n 

M n = l[m j =E(Z n \Z = l), (n>l). 

3 = 1 

The proof is (barely) indicated by Agresti (1975), and with differing notation. 

The following result is fundamental. 
Theorem 7 The conditioned process (Z n \Z n > 0) has a limit law iff linin^oo B n = B € [0, oo]. Suppose 
this condition holds, (i) If < B < oo, then 

(z n \z n >o)S~ 

where 



1 + B-Bs 
(ii) If B = oo then 

(Z n /B n \Z n > 0) 4 £, 

where £ has a standard exponential law. 
Proof Observe that 

(l + B n -A n )(l-s) 
W " (S)= (l-s)B n + l ' 

whence 
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Assertions (i) and (ii) are obvious if B exists, and the only other possibility is that (B n ) has many limit 
points, and hence there exists no limit law, with or without normalization. O 



Let T = inf{n : Z n = 0} denote the time to extinction. 
Theorem 8 (i) Q := lim^oo P(Z n = 0) = 1 iff lim^oo M n = and/or J2j>i r j/Mj = 00 • 
(ii) If Q exists and B = oo, then 

Z n /M n a A- I + (1 - I)£ 

where P(I = 1) = Q = 1 — P(I = 0), and / and £ are independent, 
(hi) If Q = 1, then P(T < oo) = 1. 



(3.2) 



Proof Since r n — R n = m n — 1, we obtain 



r 3 - R 3 



Mi 



J2(M7} 1 -M7 1 ) = 1-M- 1 & l + B n -A n = M n . 



It follows that 



P(Z n > 0) = 



Mr, 



l + B n 



M-'+^rj/Mj 

3=1 



(3.3) 



Assertion (i) follows. Assertion (ii) follows from Theorem 7(ii) and the fact that (Z n /M n ) is a positive 
martingale which thus converges a.s. with no further assumptions. Assertion (iii) follows since Z n ^> if 
Q = 1, and any a.s. convergent subsequence eventually hits zero. 



Definition 1 We say that the environments are weakly varying if M = lim^oo M n exists. If this is 
the case then we have a classification similar to the trichotomy. The environments are: 

(a) Supercritical if M = oo, i.e. if J2 n ( mn ~ -0 = °°! 

(b) Critical if < M < oo, i.e. if J2n( m n ~ 1) converges; and 

(c) Subcritical if M = 0, i.e. if J2 n ( m n ~ 1) = ~ °°- 



Keiding and Nielsen (1975, Theorem 2.2) prove (3.2) in the supercritical case. This follows from 
Theorem 8(h) since B = oo if M = oo. For example we can choose, as in the classical case, M n ~ cm™, 
where c > and m > 1 are constants and sup J<n rj — 0[(m/(l + e))"]. Clearly Q < 1, Theorem 8(h) 
holds, and these conditions can be realized. However, it is possible that Q = 1 in the supercritical case. 
This occurs in our example if p n — > sufficiently fast. Specifically, if 

p n = a/b n m n k a n — mp n 

where a > and {b n } is a positive sequence, then r n /M n ~ b n /ac. Now choose 6„ = oo. Indeed, if 
S n denotes the sum in (3.3) then P(Z n > 0) ~ S^ 1 , and we see that the probability of non-extinction 
can tend to zero arbitrarily slowly. In this case (3.2) holds with 1=1, but Theorem 7(h) shows that 
conditioning on non-extinction gives a non-degenerate limit law. 

The non-classical growth rate M n ~ n 5 is achieved by choosing a n and p n so that m n = (1 + n^ 1 ) 5 
(S > 0), and then {S n } can have a finite or infinite limit. If 5 < 1 and {p n } is non-decreasing then Q < 1 
only if p„ -> 1. 
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It is possible to have B = oo when M < oo, but in such a case (3.3) shows that Q — 1, and then 
Z n /M n -+ 0. In the critical case this occurs iff ^ n r„ = oo, that is Q = 1 iff B = oo iff r„ = oo. In 
particular r„ < oo yields a non-classical restrained growth regime in the sense that Z n A 7+ (1-7)5. 

Finally, it follows from (3.3) that Q = 1 in the subcritical case, and Theorem S(iii) applies. However 
it still is possible that B — oo, with Theorem 7(h) holding. To see this, let m < 1 in the classical scenario 
above, and observe that 

n 

B„ ~ m n ^ m i r i = X! m%r ™-i- 

3=1 i=0 

Hence, if r„ — > r < oo then B = r/(l — m). On the other hand, r„ — > oo if p„ — > and then Fatou's 
lemma shows that B = oo. 

Choosing m„ = (1 + n -1 ) -5 gives the non-classical decay M n ~ rt -<5 , and then _B„ ~ n _<5 5Z™ =1 
If r„ = n a l{n) where a is a real constant and l(x) is slowly varying, then 

'0 ifa<-l, 

B„ ~ < l(n) if a = -1, 

,oo if a > — 1. 
Hence £> = oo if — 1 < a < 0, in which case p n — > 1. 



4 Maximum family size 

Consider now the maximum n-th generation family size M n := maxi< i <^„_ 1 ^Q(n). Since M n = if 
Z n -i = 0we consider the conditional distribution function 

H n (x) := P(M n < x\Z n _ x > 0) = E [F^(x)\Z n ^\ , (4.1) 

where F n is defined at (1.2). Analogues of both Lemmas 1 and 2 can be given, but here we consider only 
the one of Lemma 2. 



Theorem 9 Suppose that lim^oc p n = p with < p < 1. 
(i) If linin^oo a n — a and lim„^oo B n = B < oo, then 

l-a(l -p)W 



lim H n (x) 



l + aB(l -p)W ' 
(ii) If limn^oo B n = oo and lini^oo log(a„B„) = a, then 

lim H n {x) = \l + e a {l-p)^ 



(x > 0). 



(x > 0). 



Proof For (i) observe that F n (x) — > 1 — a(l — p)^. The assertion follows from Theorem 7(i), (4.1), and 
the uniform convergence property of the continuity theorem for probability generating functions. For (ii), 
observe that the assertion of Theorem 7(ii) is equivalent to the limit statement 



lim E (>»-i/*»-i Zn _ x > o) = [1 + logs- 1 ] 
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and the convergence is uniform with respect to s in the interval [0 < s f < s < 1]. The assertion follows 
by setting s — s n :— F„ n ~ 1 '(x) and seeing that {s n } has the limit (2.3). O 
Our next result extends Theorem 1 to the branching process setting. The proof shows that the role 
played by v n in §2 is here played by [-B„_i] and indeed that its fractional part can be ignored. Accordingly 
we define a* = log(a„B„_i), and V denotes a random variable having the standard logistic distribution 
function L(x) = (1 — e~ x ) 1 , all real x. 

Theorem 10 Suppose that lim^oc B n = oo and for some real c 

lim p n = & lim a* n p n = 2c. (4.2) 

(i) If lim n _>oo a* = oo, then 

(pA-<|Vi>0)^V-c. 

(ii) If limn^oc a* = a (— oo < a < oo), then 

( Pn M n \Z n _ 1 >0)S(V + a)+. 

r d i 

Proof The proof is essentially the same as for Theorem 9(h), letting s n :— F„ " (x n ) where x n — 
(x + a^)/p n , just as for Theorem 1. The limit of {s n } is given by (2.7) for (i), and for (ii) by G(x) if 
x > —a, = otherwise. O 



Consideration of the proofs of Theorems 1 and 9 shows that the conditions listed in Theorem 10 are 
necessary for its conclusions. The simple approach we use for Theorems 9(h) and 10 give an obvious 
analogue for Theorem 2. 

Theorem 11 If the notation and assumptions of Theorem 2 stand with a n replaced by a* , then 
L(~f(x - 1)) = liminf P(M n — C n < x) < limsupP(A4 n — C n < x) = L("/x), (-oo < x < oo). 

n >oc n — f qq 

The various conditions in Theorems 9-11 can be satisfied, but we will show that all but one set 
is satisfied by the branching process obtained from sampling the linear birth and death process (£> t ) 
at irregular times, an example mentioned due to Keiding and Nielsen (1975). Let Z n = B tn where 
< t n < t n+ i — > too < oo- If A and \i are the birth and death rates, respectively, and d n = t n — t n _i, 
then a n = m n p n , 

X — f2 



(1 + \d n ) 1 if A = /x, 



and 



I A(M„-1) -f \ / 

{ \t n if A = fi, 

The environments are weakly varying in this case, and our criticality classification coincides with the 
standard one if too = oo, and they are critical if too < oo. The latter case is vacuous as far as Theorems 
9-11 are concerned, since B < oo and p = 1, i.e., there exists no non-degenerate conditional limit law for 
M n . So we assume that too = oo and that d n — > d < oo. 
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If the environments are subcritical, A < /i, then B < oo, p n — > p with 

< p = - ~ A < 1 & 0<m oo = e-^-^ d < 1 

(1 - ATOoo 

and a„ — > a = pm^. Hence Theorem 9(i) holds and the limit law is degenerate at the origin if d = oo, 
but not otherwise. 

If A > (i then B = oo, but Theorem 9(h) cannot hold because a = (1 — X/fi)/(l — /i/moc) > 0, whence 
a* — ► a < oo is violated. We show that the conditions of Theorems 10 and 11 can be satisfied. First, 
suppose that d = oo, in which case p = 0. In the supercritical case A > \x we further suppose that the 
second member of (4.2) is satisfied, i.e., 

2c 

lim (t n /m n ) = €[0,oo). (4.3) 

n— >oo A — fJi 

For example, this condition is satisfied with c = if d n > £(n) — ► oo where £(x) is a slowly varying 
function. If (4.3) holds then 

((M n /m n ) - (A - /i)t„| Z n _ x > 4 V - c. 

Now suppose that A = fi and t„ = n s £(n) where 5 > 1, £(a;) is slowly varying and that t„/n — > oo. 
Then (4.2) holds with c = 0, a„S„_i ~ n, and again Theorem 10(i) holds in the form 

M " ' Vi>0|AV 



A5n*-^(n)' l0g " 

On the other hand, if t n -\/t n -> r £ (0, 1) then an£?„_i — > w := r/(l — r) > 1 and a n — > a = logw. 
Hence Theorem 10(b) holds in the form 



^n-i>0) 4(V + a)+. 



(1 - r)At„ 

If A > jU, too = oo but d < oo then < p < 1 and Theorem 11 holds. 

We end by observing that analogues of results in §2 for the minimum and range can be taken into the 
branching process context, but we leave this as an exercise for the reader. 
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